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Abstract 

A novel crystalline structure of hybrid monolayer hexagonal boron nitride (BN) and graphene 
is predicted by means of the first-principles calculations. This material can be derived via boron 
or nitrogen atoms substituted by carbon atoms evenly in the graphitic BN with vacancies. The 
corresponding structure is constructed from a BN hexagonal ring linking an additional carbon 
atom. The unit cell is composed of 7 atoms, 3 of which are boron atoms, 3 are nitrogen atoms, 
and one is carbon atom. It behaves a similar space structure as graphene, which is thus coined 
as 5-B3N3C. Two stable topological types associated with the carbon bonds formation, i.e., C-N 
or C-B bonds, are identified. Interestingly, distinct ground states of each type, depending on C-N 
or C-B bonds, and electronic band gap as well as magnetic properties within this material have 
been studied systematically. Our work demonstrates practical and efficient access to electronic 
properties of two-dimensional nanostructures providing an approach to tackling open fundamental 
questions in bandgap-engineered devices and spintronics. 



PACS numbers: 73.22.Pr, 75.70.Ak, 78.67.Wj 



I. INTRODUCTION 



Two-dimensional (2D) nanomaterials, such as graphene and monolayer hexagonal boron 
nitride (h-BN), are expected to play a key role in future nanotechnology as well as to pro- 
vide potential applications in next-generation electronics. Recently, novel hybrid structures 
consisting of a patchwork of BN and C nanodomains (BNC) was synthesized through the 
use of a thermal catalytic chemical vapour deposition method 1 . This finding immediately 
has attracted a great deal of research interest^-, given it demonstrating a hitherto efficient 
route to tune the band gaps of these 2D materials. 

It is well known that the perfect hexagonal and planar structure of BNC largely depends 
on the well matching between BN and C domains. However, it is indeed an outstanding 
challenge as BN and C phases are naturally immiscible in 2D 1 . This explains why Ci et a/.- 
could observe some wrinkles in atomic force microscopy (AFM) image. Such mutual contra- 
diction mainly originates from the domain boundary effect and the staggered potentials of 
B and N atoms in BNC, which doubtlessly affects their continuous tunable electronic energy 
gaps. It has been confirmed in the related theoretical calculations^- where the band gaps 
behave a strong oscillation feature. 

The X-ray photoelectron spectroscopy (XPS) of BNC measured in the work from Ci et 
a/.- shows two additional types of C bonding configurations, which correspond to the C-B 
and C-N bonds with the bonding energies of around 188.4 and 398.1 eV, respectively. This 
feature means that two inequivalent C bonding types, i.e., C-B or C-N bonds, must be 
present in the boundaries of the hybridized BN and C domains, which have a significant 
effect on determining the local structures and subsequently vary the electronic properties in 
this system. For example, the transport channels behave a robust characteristic gap when 
the topological index changes sign of the valley Hall effect^. In addition, both Raman D 
band at 1,360 cm -1 and D' band at 1,620 cm -1 are also observed in BNC 1 , which were 
attributed to the lattice disorder or the finite crystal size. This lattice disorder effect might 
directly introduce vacancies to these 2D hexagonal system^— , which is also true in BNC 
as shown in Fig. 2a of Ref.-. 

For the planar BNC structure, although the larger domains would be preferred to de- 
crease the total domain interfacial energy, the randomly distributed hybrid domains, the 
immiscible phases as well as the induced vacancies must result in various complex structures 
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FIG. 1. The total energy per atom as a function of lattice constant a for (7-B3N3C. Yellow, green, 
and gray balls represent B, N, and C atoms, respectively. Their respective 2x2 supercells are also 
given nearby. 

of BNG^~— . As an immediate consequence is the largely inaccessible synthesis of expected 
BNC in experiments 17 , which next hinders realization of bandgap-engineered applications in 
actual devices. In this endeavour, exploring the structures and electronic properties associ- 
ated with C bonds formation in BNC that contain C-N or C-B bonds is an interesting topic 
that must be addressed before widespread synthetic applications. Thus, a simple model 
of BNC where the C bonds play a crucial role must be considered again. More impor- 
tantly, a deep theoretical understanding, which originally was concealed behind the complex 
hybridized structures, is imperative. 

Here we report such a simple model of BNC may be just the graphitic B3N3C ((7-B3N3C), 
which is a perfect 2D monolayer graphite-like structure as shown in Fig. 1. As mentioned 
above, C atom can bond to B and N atoms. Therefore, two topological types of g-B 3 N 3 C 
are easily deduced. One is the a-g-B 3 N 3 C related to the the higher bonding energy of 
C-N bond, while the other one is the ft-g-^^N^C constructed based on the lower bonding 
energy of C-B bond. It can be seen that such a material is essentially C-doped graphitic 
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BN (g-BN) with vacancies*^. The substitution of the N/B atom with a C atom in g-BN 
with the B/N vacancy will yield the a//3-g-B 3 N 3 C structure. The interactions among the C 
atoms and/or vacancies as well as the C bonding types (C-B and C-N bonds) in g-B 3 N 3 C 
significantly alter its electronic properties. To explore this effect, standard density functional 
theory with different functionals (see the following discussion) calculations have been carried 
out for this predicted material. Remarkably, such material displays two distinct electronic 
structure properties: a-g-B 3 N 3 C is a semiconductor, while /3-g-B 3 N 3 C behaves metallic and 
leads to a magnetic ground state. 

This paper is arranged as follows: In the second section, we present the computational 
method used in this work, followed by electronic band gap in a-g-B 3 N 3 C and magnetism in 
/3-g-B 3 N 3 C in the third section. We then conclude this paper in the fourth section. 

II. COMPUTATIONAL METHOD 

For structural optimization, we employed density functional theory with the general- 
ized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) 1 - for the exchange- 
correlation (XC) potential within the projector augmented wave method as implemented 
in VASP 19 . An all-electron description, the projector augmented wave method, is used to 
describe the electron-ion interaction. The cutoff energy for plane waves is set to be 500 
eV, and the vacuum space is at least 15 A, which is large enough to avoid the interaction 
between periodical images. A 7x7x1 Monkhorst-Pack grid is used for the sampling of the 
Brillouin zone during geometry optimization. All the atoms in the unit cell were allowed 
to relax, and the convergence of force is set to 0.01 eV/A. Additionally, spin polarization is 
turned on during the relaxation processes. All other calculations of accurate electronic prop- 
erties were performed using the full-potential linearized augmented plane- wave (FLAPW) 20 
method as implemented in the WIEN2k code 21 . It is well known that different XC potentials 
can lead, depending on the studied materials and properties, to results which are in very bad 
agreement with experiment, e.g., for the band gap of semiconductors and insulators which is 
severely underestimated or even absent^ 2 -. For this reason, the modified Becke and Johnson 
(MBJ)22i2i potential in the framework of local-density approximation (LDA}2£>2& is taken to 
calculate the band gap of a-g-B 3 N 3 C, while the magnetism in /3-g-B 3 N 3 C is described by 
the PBE (the standard GGA for materials) potential. 
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III. RESULTS AND DISCUSSION 



A. Electronic band gap in a-g-B 3 N 3 C 

Figure 1 plots the total energy per atom against the lattice constant a (the lattice constant 
c is fixed) for g-B 3 N 3 C. It can be seen that the total energy of g-B 3 N 3 C as a function of lattice 
constant a has a single minimum, meaning that the geometrical structure would be stable. 
Particularly, the charge population analysis reveals that the electron density around the C-N 
bond in a-g-B 3 N 3 C is much higher than that on the C-B bond in /3-g-B 3 N 3 C, showing that 
the C-N bond is relatively strong, which is also consistent with the experimental results 1 . 
This strong interaction between C and N atoms in a-g-B 3 N 3 C directly results in a short 
C-N bond length (see following), and can balance the strain of the monolayer graphite-like 
structure. Thereby a-g-B 3 N 3 C could be a more thermodynamically stable topological phase 
against /3-g-B 3 N 3 C, which has a lower total energy of around 0.43 eV/atom compared with 
^-B 3 N 3 C. 

To explore further the mechanical stability of a-g-B 3 N 3 C, the optimized lattice constant 
a = 4.67 A is first obtained, as depicted in the left panel of Fig. 1. Importantly, the 
C-N bond length converged to 1.31 A, considerably reduced from the typical bond lengths 
1.37-1.48 A in the related materials^. This strongly suggests that the nature of the higher 
binding energy of C-N bond^, which also influences the B-N bonding and extends its length. 
The obtained B-N bond length is 1.48 A, which is slightly bigger than the value of 1.45 A 
in h-BN. The calculated partial density of states (PDOS) is shown in Fig. 2. The valence 
band is dominated by B p and C p states, while the conduction band is only dominated by 
B p states. There one can find that the majority of C p states to be semicore lying 6-9 eV 
below the Fermi level. These states interact with those comprising the valence band with 
the same symmetry. As a result there is a small admixture of C p and N p states close to the 
Fermi level. However, in a-g-B 3 N 3 C one finds significant admixture of C p and N p states 
in the semicore energy window with 6-9 eV below the Fermi level. This suggests that the C 
p-N p interaction in semicore region contribution to the C-N bonding is significantly more 
important in a-g-B 3 N 3 C than the C p-N p interaction close to the Fermi level. 

Now, let us look at the band structures of a-g-B 3 N 3 C, as given in Fig. 3. It explicitly 
demonstrates that all three XC potentials give similar band structures. The highest occupied 
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FIG. 2. Total and partial DOS of a- 5 -B 3 N 3 C for three XC potentials LDA, PBE, and MBJ. The 
vertical dotted line denotes the Fermi level, and also indicates end of the fundamental band gap 
which starts at E — E-p = eV. The black, blue, and red lines correspond to the total, p, and s 
DOS, respectively. 



crystalline orbitals (HOCOs) are located at the K point of the reciprocal space, while the 
lowest unoccupied crystalline orbitals (LUCOs) appear at the M point. This leads to an 
indirect-band-gap semiconductor. To obtain the band gap more accurately [see Fig. 3(c)], 
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FIG. 3. Calculated band structures for a-g-H^N^C model with the three XC potentials LDA (a), 
PBE (b), and MB J (c). 

the MBJ XC functional is used^ 1 ^. The band gap of a-g-B^N^C within MBJ is obtained 
to be 1.22 eV, which nicely locates the middle region between 0.59 and 1.80 eV for the 
BNC samples with 12.5% and 50% C contents, respectively^. Note that the C content 
in a-g-B 3 N 3 C is 25%. In the experiments^, the absorption edges are redshifted as the C 
concentrations increase, which shows a tunable mechanism of optical band gap in actual 
applications. By comparing with Ref.- where BNC with around 65% C concentration shows 
an optical bandgap of 1.62 eV, we infer that such a higher energy absorption edge (take 
into account the band gap of 1.22 eV in our case with 25% C concentration) arises from 
the formation of individual BN and graphene domains. In this way, the even distribution 
of C in BNC systems might serve as a good guide to find alternative solutions to existing 
bandgap-engineered applications. Future researches can test this prediction directly. 

In addition, the band gaps based on the LDA and PBE methods are equal to 0.83 eV 
[see Figs. 3(a) and (b)]. The relative band gap correction from MBJ with respect to 
LDA, TZ = (A MBJ - A LDA )/A MB j with A MBJ / LDA being the band gap, is about 32%. We 
can see that our calculated 71 for a-g-B 3 N 3 C ( shows an excellent agreement with the 16 sp 
semiconductors) lies within the range 16.0%-100.0% (see Table 1 for details). The correction 
value is particularly very close to that in BN (25%), GaN (42%), A1P (37%), and A1N (25%). 
This finding is not surprised because the listed 4 solids have at least one element close to that 
of a-g-B 3 N 3 C in the periodic table of elements. It means the similar chemical circumstances 
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TABLE I. The band gaps (in eV) and the relative band gap correction 7Z (%) of 16 sp semicon- 
ductors and the predicted a-g-BsNaC. The theoretical and experimental band gaps of the 16 sp 
semiconductors are directly taken in2^. 
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in these materials can be well described by the same XC functionals. This underlyingly 
confirms our prediction validity and the a-g-B^N^C might be carried out experimentally. 

From Fig. 2, which shows the DOS of a-g-B^N^C, we can see that the effect of the MBJ 
potential is to shift up (with respect to LDA/PBE) the unoccupied B 2s and 2p states. 
Here, three major differences between the LDA/PBE and MBJ methods can be extracted, 
(a) Another obvious effect of MBJ potentials is to shift down the middle of the valence 
band at around -4.0 eV. (b) The hybridization of s and p states of dominant B 2s and 
other atoms 2p states at the bottom of the valence band is very strong in MBJ calculations 
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FIG. 4. The optical absorption expressed by the imaginary part of the dielectric tensor e<i for 
a-<7-B 3 N 3 C. The imaginary part of the dielectric tensor is averaged for four different directions. 
Three XC potentials including LDA (a), PBE (b), and MB J (c) are shown in calculations. 

(denoted as down arrows in Fig. 2). (c) The correction character of a-g-B 3 N 3 C is much 
more pronounced with the MB J than with the LDA/PBE, which narrows the valence band 
just below the Fermi level. We would like to stress that the MBJ potentials open a band 
gap of 0.39 eV in the a-g-B 3 N 3 C model compared to the result of LDA, which is consistent 
with the orbital-dependent potentials principle^. 

The optical absorption spectrum for a-g-B 3 N 3 C within MBJ potentials [see Fig. 4(c)] 
shows a big absorption packet with two adjacent peaks in the range of 2.5-5.0 eV, which 
originates from the band structure as shown in Fig. 3(c). As an indirect band gap material, 
the transition from the K to T point is very weak as the momentum conservation rule 
is not satisfied here, and thus the corresponding characteristic photoluminescence in the 
optical absorption spectrum can be negligible. This is also true in our case a-g-B 3 N 3 C. The 
first peak corresponds to the direct band gap transition at the M point. The second peak 
comes from the larger direct gap from higher energy states located at the T point. More 
importantly, these two peaks are not separated distinctly. This feature can be attributed 
to the even distribution of C atoms in a-g-B 3 N 3 C as mentioned above, which behaves a 
different formation mechanism compared with the hybrid BNC in experiment^. 
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Correction effects were taken into account by adding a LDA correction potential in MB J2^. 
This important physical effect opens an additional band gap through mimicing very well 
the behavior of orbital-dependent potentials and causes a rigid blueshift of the absorption 
spectrum compared with the LDA/PBE curves as shown in Fig. 4. This explains the 
excellent qualitative agreement of the hybrid exchange potential optical absorption spectrum 
seen in Fig. 4, due to a compensation of significant errors within the standard DFT methods. 

B. Magnetism in /3-#-B 3 N 3 C 

The predicted structure of /3-g-B 3 N 3 C is shown in the right panel of Fig. 1. The lattice 
constant of /3-g-B 3 N 3 C is obtained to be 5.26 A as shown in Fig. 1. The results show that 
the equilibrium value g?bc = 1-52 A, which is close to the value of graphite-like BC3, namely, 
1.55 A—. It is to be noticed that all the equilibrium values cIbn in a-g-B 3 N 3 C are equal to 
1.42 A, which is slightly less than the value of 1.45 A in pristine BN sheet. This implies the 
stronger B-N bonds formed in /3-g-B 3 N 3 C Our calculations show that the /3-g-B 3 N 3 C leads 
to a ground state with a magnetic moment of 0.68 /ib- The nonmagnetic state is 0.07 eV 
higher than this ground state. 

From the band structures we see that although both the pristine BN sheet and graphene 
are nonmagnetic, the /3-g-B 3 N 3 C model can be spin polarized. It is necessary to discuss 
magnetism in more detail from its electronic structures. Figures 5(a)-(c) present the band 
structure and spin-resolved total density of state (TDOS). Remarkably, two bands cross 
the Fermi level [black and red circles in Fig. 5(a)], and make the Fermi energy level oc- 
cupied completely along the entire high-symmetry lines. In contrast, the spin-down one 
does not possess such a strongly localization feature, just one band (black circles) crosses 
the Fermi level monolithically as shown in Fig. 5(c). The calculated magnetic moment in 
/3-g-B 3 N 3 C should originate from this asymmetric spin-dependent localization. The corre- 
sponding strong spin splitting can be further confirmed from TDOS as shown in Fig. 5(b). 
In addition, close examination of the top valence bands [see Fig. 6(a)] indicates that the 
strong localization mainly comes from the 2p atomic orbitals of C and N atoms. Here we also 
show that the Fermi surfaces of /3-g-B 3 N 3 C in the first Brillouin zone. The unique feature 
of Fermi surface almost parallel to the high-symmetry lines [see Figs. 5(d)-(f)] is the direct 
manifestation of the bands near the Fermi level in Figs. 5(a) and 5(c). 

10 




\ t 










r m k 
Band structure 



DOS 



M K r 
Band structure 




FIG. 5. (a)-(c) Band structures and spin-resolved TDOS for /3-(7-B 3 N 3 C. The dotted line indicates 
the Fermi level. The arrow denotes the spin polarization direction: up for spin up and down 
for spin down, (d)-(f) Fermi surfaces drawn in the first Brillouin zone and the corresponding 
high-symmetry points, (d) and (e) for spin up and (f) for spin down. 



It should be notice that such magnetism is induced without transition metals and with- 
out external perturbations, so that /3-g-B 3 N 3 C behaves as the first, theoretically predicted, 
metal-free magnetic material in hybrid BNC system. Apparently, our finding points out a 
new direction for further related experimental investigations in spintronics. 

We now address the possibility to magnetism origin in /3-g-B 3 N 3 C. Based on the orbital- 
resolved density of state as shown in Fig. 6(a), the magnetic moment are mainly ascribed to 
the 2p orbitals of C and three N atoms. The spin polarization of C atom offers a magnetic 
moment of 0.18 /xb- 0.18 /Jb of total magnetic moment are shared equally by the 2p orbitals 
of three N atoms. The remaining magnetic moments distribute evenly in the interstitial 
region among N atoms. This is conceivable, because /3-g-B 3 N 3 C, if compared with other 
BNC systems, has large interspaces between N atoms [see the 2x2 supercell in Fig. 1 or 
Fig. 6(b)]. The remarkable feature is that the equilibrium surface density of /3-g-B 3 N 3 C is 
around 1.27 times larger than a-g-B 3 N 3 C. Owing to the large interspaces between atoms, 
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FIG. 6. (a) PDOS on the interstitial region and on the 2s and 2p orbitals of C and N atoms in 
/3-g-B 3 N 3 C. (b) The 3D iso-surface plot of spin density for the 2x2 supercell at the value of 0.04 
e/A 3 : (c) A side view of spin density corresponding to (b). 



the hydrogen storage in /3-g-B 3 N 3 C may be expected^. The detailed description of hydrogen 
storage related to /3-g-B 3 N 3 C is beyond the scope of this work. 

Concerning the special carbon atom in /3-g-B 3 N 3 C, the counting of four valence electrons 
is as follows: Three electrons participate in the sp 2 hybrid orbital, which forms a planar 
structure. The remaining one electron is then redistributed in the whole unit cell due to the 
enhanced B-N covalent bond (with shorter bond length compared with the value in pristine 
h-BN), which makes the magnetic properties more complicated. From the fourth electron, 
only 18 percent of the electron still fill the 7r-orbital of C atom and contribute a magnetic 
moment of 0.18 /xb- This is in fairly good agreement with theoretical description of the ir- 
orbital state to the local magnetic moment ~0.3 /ib in the graphene-like systems 3 ^ 1 ^. Around 
32% of the fourth electron mainly reside at the interstitial region, which play a crucial role 
in /3-(yf-B 3 N 3 C as: (a) Enhance the B-N covalent bond, (b) Provide the main interstitial 
magnetic moment of 0.32 /ib- (c) Promote the 2p z of N atom spin polarized slightly with 
a magnetic moment of 0.06 /iB per N atom. The remaining of the fourth electron act as 
the conduction electrons and make the system metallic, which dominate the mechanism of 
ferromagnetic ordering in /3-g-B 3 N 3 C. In the case of /3-g-B 3 N 3 C, we can see that the electron 
spin at the localized 7r-orbital state of C and N atoms, as well as the interstitial region compel 
two energy bands localized strictly along the entire high-symmetry lines, i.e., the T-M-K-r 
line [see Fig. 5(a)]. Thus, the RKKY interaction 3 ^ 3 - among the magnetic sites through the 
residual conduction electrons forms a spin ordering in these orbitals, which is the physical 
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origin of ferromagnetism in P-g-B^N^C Figures 6(b) and 6(c) respectively plot the top and 
side views of the 3D iso-surfaces for net magnetic charge density in xy plane. This finding 
is insightful and three major points deserve comment: (a) The C site is more spin polarized 
as compared with each N site, (b) The dumbbell-like magnetic moment distribution along z 
direction implies the 2p z orbital becomes partially filled with one spin-up electron, (c) The 
induced moments are ferromagnetic coupled between the N and C sites based on the RKKY 
exchange interaction model as mentioned above. 

IV. CONCLUSIONS 

In summary, we have predicted a novel crystalline material g-B 3 N 3 C, which displays two 
distinct electronic properties where the selective bonding type of C atom is a key parameter 
for future industrial processes. a-g-B^N^C is a semiconductor, while jS-g-B^N^C behaves 
metallic and holds a magnetic moment of 0.68 /zb- Importantly, compared with the hybrid 
BNC, (7-B3N3C is proposed to have a simple structure, which can be applied in various fields 
due to its unique properties. 
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